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Abstract 



< 

5^ In this note we analyze a model for a unidirectional unsteady flow of a 

ri viscous incompressible fluid with time dependent viscosity. A possible way 

I " I to take into account such behaviour is to introduce a memory formalism, 

including thus the time dependent viscosity by using an integro-differential 

^-H term and therefore generalizing the classical equation of a Newtonian 

^ viscous fluid. A possible useful choice, in this framework, is to use a 

r ■ rheology based on stress/strain relation generalized by fractional calculus 

modelling. This is a model that can be used in applied problems, taking into 
account a power law time variability of the viscosity coefficient. We find 
. analytic solutions of initial value problems in an unbounded and bounded 

^^ domain. Furthermore, we discuss the explicit solution in a meaningful 

^^ particular case. 
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"k^ 1 Introduction 

5—1 

C^ In recent years fractional calculus modelling has found many interesting ap- 

plications in the description of viscoelasticity (see for example Mainardi [T]). 
It is well known that replacing the ordinary time derivative with the Caputo 
fractional derivative in the constitutive rheological equations we introduce a 
memory formalism in our problem. The Caputo fractional derivative is defined 
by a convolution of the ordinary integer-order derivative with a power law kernel. 
In this sense, it is evident the utility of a generalised rheology for complex 
media. Within this framework, several applications of fractional modelling of 
non-Newtonian fluids were recently investigated in many papers. For example 
Tong et al. |2j treated the unsteady helical flow of a generalised Oldroyd-B 



fluid and Qi and Xu [3] studied the Stokes's first problem for such generalised 
Oldroyd-B fluid. Moreover, Vieru et al. |4j studied the fractional Maxwell model 
and Hyder Ali Muttaqi Shah j^ studied the fractional Burgers's fluid model. 

In this paper we discuss a model of unsteady unidirectional flow of an 
incompressible fluid with time-dependent viscosity in a simple geometry. The 
main idea is that a possible way to take into account such variability is to 
generalise the equations of a viscous Newtonian fluid by introducing an integro- 
differential term, i.e. a memory formalism. Using this approach we include the 
time-dependent viscosity in the kernel of this term thus suggesting an analytic 
way to understand the physical role played by the variable viscosity. 

A possible useful choice is to use a rheology based on a stress/strain relation 
generalised by means of fractional calculus. We arrive then at a class of integro- 
differential equations based on a power law kernel. This analytic tool provides 
us useful mathematical instruments that can be used in applied problems which 
consider a viscosity coefficient with a power law time variability. 

We show that this fractional model of an unsteady flow is described by an 
integro-differential equation with a Riemann-Liouville fractional integral. We 
also discuss a general way to solve the initial value problem related to the master 
equation in an unbounded and bounded domain. Finally, we discuss the explicit 
solution in an interesting case in order to better understand the physical meaning 
of our results. 

2 Unsteady unidirectional flow of incompress- 
ible fluids with time- dependent viscosity by 
fractional calculus modelling 

The classical stress/strain relations of a Newtonian viscous fluids are defined by: 

Tij = 2^eij + ATr(eij) (2.1) 

= 2^e,j+A(V-F), 

where tij is the Cauchy stress tensor, /x is dynamic viscosity of the medium, e^ 
is the strain tensor and A is the volume viscosity coefficient. As we are studying 
incompressible fluids, that is when \/ ■ V = 0, then the second term in the right 



hand side of (2.1 1 vanishes and we have the simple relation 

Tij = 2fieij. (2.2) 

In this note we suggest a way to take into account a possible time-dependent 



viscosity by generalising (2.2 1 with an integro-differential relation: 



, f K{t-T)e.,j{T)dT, t>0 
Jo 



r„ = 2^0 / K{t - T)e^jiT) dr, t > 0, (2.3) 

Jo 

where fiQ is the constant viscosity of the specific memory less case and K{-) is 
a general kernel; its form depends on the type of time variability of the fluid 



viscosity. We notice that equation (2.31 was first introduced by Gerasimov [5] 



in relation to the application of generalised laws of deformation to problems of 
internal friction (see also Mainardi |Tj and Uchaikin ^Tj, pages 285-286). 



This proposal has a clear physical meaning: the role played by the viscosity 
time dependence is to introduce a memory formalism and hereditary effects on 
the fluid flow. 



We now propose a useful way to generalise (2.2), specialising (2.31 with 
a power law kernel. This generalisation is based on the application of the 
Riemann-Liouville fractional integral. Hence, we have 

r„- = 2^0 Jre»j, a e (0, 1), t > 0, (2.4) 

where J" is the Riemann-Liouville time fractional integral defined by [8] 

Jtfit) = ^^j{t-Tf'^f{T)AT, /3>0, t>0. (2.5) 



Within our framework, the meaning of (2.4) is the following: through fractional 



integration we take into account a power law time variability of the viscosity 
coefficient i.e. we include non-local effects. This can be interesting in relation 
to the study in the field of time variable viscosity in general (for example in 
solid suspensions lOj). Here we recall that there are many different approaches 
to fractional rheologies starting from the seminal works of Scott-Blair and 
Gerasimov (for a historical review the reader can refer to Mainardi PD]). However 
we are not discussing the rheology of a specific class of viscoelastic media, we 
used the fractional formalism only in order to include the time variability of 
the viscosity by means of convolution. Therefore our model cannot be directly 
related to the already classical fractional rheologies widely discussed in the 
literature (see for a recent review Surguladze [H]). Instead we are going to study 
the ffow of a viscous fluid with a power-law time dependent viscosity. In general 
it leads to an integro-differential generalisation of the Navier-Stokes equation. 
In the following we will discuss a particular case as the more general case leads 
to a complex nonlinear integro-differential equation which however could be the 
object of a further investigation. 

The equations of flow of an incompressible fluid without volume force in 
Cartesian coordinates are 

where D/Dt = 9f + (F • V) is the material derivative. For simplicity, we assume 
a constant density p — 1- We treat the unidirectional motion of this fluid on 
a plate. The flow of a unidirectional velocity field, in Cartesian coordinates is 
given by 

V^u{v,t)i, (2.7) 

where u is the velocity in the cc-direction, z is the unit vector in the x-direction, 
X is the coordinate along the plate and y is the coordinate normal to the plate. 
Considering these hypotheses on the velocity field, taking Vp = (constant 



pressure), and using (2.4) we have that 

dtu = ^iQdy{J^dyu), t>0. (2.8) 

This is an integro-differential equation with a fractional Riemann-Liouville 
integral. 



In the next, we first discuss the analytic solution of the following initial value 
problem (IVP) in an unbounded domain. The following step will be, in Section 
|3] the more realistic case of a bounded domain. 



dtu = fiodyiJj^dyu), y e K, t>0, 



(2.9) 



In order to solve this IVP we make use of the Laplace-Fourier transform. 
Recalling that the Laplace tranform of the fractional integral is given by [5] 



C{jru{y,t)){s) = s-"u{y,s), 



(2.10) 



where s is the Laplace transform parameter and u[y, s) is the Laplace transform 
of u(jj,t), we immediately obtain 



su{uj, s) — f{i^) = —i^ Mos "u(a;, s), 
where u{(^, s) is the Laplace-Fourier tranform of u{y, t). Therefore 

s"/(w) 



u{uj, s) 



«Q + 1 



A*o^ 



We now recall the following useful relation [5]: 



e~'H''-^E^^^{at'')dt = 



where 



E,^,{z)^Y. 



— ' Tii>r + li) 



z e 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



is the Mittag-Leffler function. The Laplace transform in (2.12) can be easily 
inverted now, leading to the Fourier transform u{uj,t): 



uiLj,t) = f{io)E^+i^i{-fio^h"+'). 



(2.15) 



By recurring to relations (14.7) and (14.8) of Mathai and Haubold J^ we can 
invert also the Fourier transform. With this inversion we can write the solution 
u{y,t) by means of a Fox function. Indeed we have 



u{y,t)^f{y)^ 'M^1'° 



\y\ 



Hot° 



(l,a + l) 
(1,2) 



i > 0, J/ e K, (2.16) 



where * is the Fourier convolution. By means of simple manipulations and 
considering the representation of the Fox function as a contour integral [12j we 



obtain that 


u{y,t) 


= /(y) * 




= /(y) * 




= /(y) * 



^^ f r(l + 2z) / \y\^ 



dz 



(2.17) 



1 



1 



r(i + 2z) 



-2z-l 



\y\ 



^1/2^^ 27rz 7c r(l + {a + l)z) {^^l^h^i^ J 
1 1 f r(r) / |y| 



dz 



/(y) 



2^.5' ^t- 
1 



i/2^-±i 2^z 7c r (1 + (a + 1 



'^; 






1 



r(r) 



|y| 



r, 1/2, ^Hzi 0~j / p /l Q + 1 I Q + l_\ I 1/2,5±1 



dr 



dr 



Wright special function, defined also as 



In the last step of (2.17) we recognize the Mellin-Barnes representation of the 

(2.18) 



W,,^{z)^Yl 



r^ r-!r(i/r + fi) 



v > —1, fi e 



z e 



We therefore arrive at the final following form. 
1 



uiy,t) = /(y) * 



„ 1/2, ^H;^ 



-w 



m 






t>0, y eR. 
(2.19) 



Note that the above solution coincides with the solution to the superdiffusive 
fractional diffusion equation with the initial condition dtu{y,t)\t=o = 0. Note 
also that the well-known fundamental solution to the fractional diffusion equation 
written in terms of Wright functions (see e.g. Mainardi [13 ) can be retrieved by 
letting f{y) = S{y). 

This result is apparently not surprising, but it forces us to make some 
considerations on a sensitive topic related to fractional calculus, i.e. the validity 
of the semigroup property for fractional derivatives. Actually, going back to 



the master equation (2.8) and applying the Caputo fractional derivative to both 
sides, we obtain 

d^dtu = fiodyyU, t > 0, y e M. (2.20) 

We recall that the definition of the Caputo partial fractional derivative 9f , 
ae (0,1), is E] 



drgiy,t)^Jl-'^dtgiy,t) 



(2.21) 



r(l-a) 



{t-Trdrg{y,T)dT. 



Equation (2.20) is a sequential fractional differential equation, different from the 



fractional superdiffusive equation thoroughly treated in literature. The reason of 
the difference is that the semigroup property is not always satisfied by Caputo 
fractional derivatives, i.e. d^d^u ^ ff^^^u ^ d^d^u. 



In literature, several attempts to solve equations with fractional sequential 



derivatives as in (2.20) are present (see for example Podlubny [S]), but they 
require initial conditions with a not clear physical meaning. However, in our 
case we have a well posed problem because the master equation is basically a 
first order integro-differential equation. The difference with the superdiffusive 



fractional case is in the initial conditions. In order to solve ( 2.8 1 it is sufficient one 



initial condition, while in the superdiffusive case two conditions are mandatory. 



We remark, however, that the analytic solution of (2.8) coincides with that of 
the superdiffusive fractional equation in the specific case dtu{y,t)\t=o = 0. 

In the next section we outline a general way to solve analytically the master 
equation in a bounded domain. Moreover, we discuss the exact solution for an 
interesting application which allows us to highlight the physical meaning for this 
model. 



3 Exact solution to the master equation in a 
bounded domain 

A number of papers are devoted to the analytic investigation of exact solutions 
to fractional integro-differential equations in bounded domains. Here we recall 
for example that Agrawal [14; found a general way to solve the fractional 
diffusion equation in a bounded domain by using Laplace and finite sine tranform 
techniques. 

Here we are going to make use of similar methods to solve the following 
initial/boundary value problem: 

dfU == fiaJtdyyU, a e (0, 1), 

u(0,f) =-u(/i,i) =0, t > 0, (3.1) 

"(y,0) = /(y), 0<y<h. 



The finite sine transform of (3.1 ) is simply given by 

dtu{n, t) = /ioa^n^ Jf u(n, t), (3.2) 

where n is a wave-number and a = ir/h. Note also that u{n,t) is the finite sine 
transform of u{y, t) and is defined as 

u{n,t) ~ / u{y , t) sm{a n y) dy . (3-3) 

Jq 

Therefore, the related initial value problem (in wave-number domain) is 

\dtu{n,t) = ^oa'^n'^J^u{n,t), 
\u{n,0) = J^ f{y)sm{any)dy. 



By taking the Laplace transform of (3.4), we obtain 

~ s"u(n, 0) 



"+i+/io(H'' ^^'^^ 




10 
Time (arb. units) 



Figure 1: The Mittag-Leffler function E^^i{~f^) for 7 = 1.2 (black), 7 = 1.5 
(red) and 7 = 1.8 (green). 



Finally we can write that 
2 °° 



a^n^r+i) 



sin(arij/) / /(z) sin(anz)dz, (3.6) 




where i > 0, y G [0, /i]. 

This is a general way to solve the fractional diffusion equation in a bounded 
domain, with the additional constraint of homogeneous boundary conditions. 



Now we discuss an application of this analytic way to solve (3.1 1 in a simple 
case. We study the following problem: 



dfU — fiodyyJ^u, a e (0, 1), 

■u(0,i) =u(/i,i) =0, t>0, 
u(2/, 0) = Uo = const, < y < h. 



(3.7) 



The physical meaning of these conditions is clear: we have an initial constant 
velocity field u{y, 0) = Uq in a confined plane, with null boundary conditions at 



both sides. By using (|3.6|), we find 
u{y,t) 



Uo 



2^ 



n=l 



„ / , , o „, IS sin(anw) , , s, 

Ea+i (-V^to a^nh"+^) — !^ ^ (1 - cos(n7r)) 



(3.8) 



We have to notice that the Mittag-Leffler function of order 1 < (a + 1) < 2 is also 
called in literature as a fractional oscillation (see Fig. [I]) . The result of memory 
effects due to time variability of the viscosity coefficient is thus an anomalous 
fractional oscillation. Furthermore, we can recover the exact periodic solution 
for a = 1. As expected, this is the solution that can be found by separation of 
variables from the ordinary wave equation. 

Moreover, this solution suggests that fractional modelling is a good instrument 
to describe global damping due to local time variability of the viscosity coefficient. 



4 Conclusions 

In this note we suggested the use of a generahsed stress/strain relation to consider 
a time variabiUty of the fluid viscosity coefficient. Treating a unidimensional 
model, we showed that, under reasonable mathematical hypotheses, the fluid 
flow is described by a linear integro-differential equation involving a Riemann- 
Liouville fractional integral. We also discussed, by means of an exactly solvable 
example, the meaning and utility of this model also in relation to the study of 
the behaviour of complex fluids with time variable viscosity. Moreover, we notice 
that the real order of integration is a free parameter that can be exploited in 
order to have a better agreement with the experimental data. It is also relevant 
to understand the physical meaning of this real order of integration in relation 
to the characteristics of the fluid. Finally, we remark that within this framework 
we take into account possible emerging empirical power law behaviours by means 
of simple analytic models. 
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